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Abstract 

We study numerically the deformation of sessile dielectric drops 
immersed in a second fluid when submitted to the optical radiation 
pressure of a continuous Gaussian laser wave. Both drop stretching 
and drop squeezing are investigated at steady state where capillary 
effects balance the optical radiation pressure. A boundary integral 
method is implemented to solve the axisymmetric Stokes flow in the 
two fluids. In the stretching case, we find that the drop shape goes from 
prolate to near-conical for increasing optical radiation pressure what- 
ever the drop to beam radius ratio and the refractive index contrast 
between the two fluids. The semi-angle of the cone at equilibrium de- 
creases with the drop to beam radius ratio and is weakly influenced by 
the index contrast. Above a threshold value of the radiation pressure, 
these "optical cones" become unstable and a disruption is observed. 
Conversely, when optically squeezed, the drop shifts from an oblate to 
a concave shape leading to the formation of a stable "optical torus". 
These findings extend the electrohydrodynamics approach of drop de- 
formation to the much less investigated "optical domain" and reveal 
the openings offered by laser waves to actively manipulate droplets at 
the micrometer scale. 
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1 Introduction 



Since the seminal works of Zeleny (IZelenvl. 1 1 9 



j) an d Taylor (jTavlorl . [l964) 



1895) on dielectri c liquid rise 



on electrified jets and drops, of Pellat (jPellatl . _ 

betw een vertically oriented parallel electrodes, or of Lippmann ( Lippmaml 
mi) on electrowetting, the control of fluid interfaces deformation under an 



applied field has received an increasing attention for the past century. Such a 
control technique plays a significant role in emergi ng micro / nano t echno logies 
with applica tions as different as electros praving (jGarian- Ca.lvcl . l2004h . ink- 
j et printing (Baddie and de Lange . 19971 ). electrospin ning of polymer fibre s 



(Hohman et al 



2001al [bj), surface relief patterning (ISchaffer et all 2000) 



tunea ble lensing ( Berge and Peseux , 2000l ) , or microfluidics ( Squires and Quake , 
20051 ). When acting on finite volume objects such as drops, the electric stress 



tends to elongate the fluid interface in the direction of th e electric field, as de- 
scrib e d by th e well- e stablished ele c trohydrostatic theory JGarton and Krasucki . 



1964 : iTavlorl . Il964 : iMiksisl . Il98ll : iBasaran and Scrivenl . ll994al Jbh. This ap- 



proach is usually relevant to the extreme situations where either both fluids, 
inside and outside the drop, are insulating dielectrics with no free charges 
present at their interface or the drop is a conducting fluid while the sur- 
rounding one is insulating. For a drop with zero net electric charge placed in 
a uniform external electric field, the electrohydrostatic theory predicts that 
the drop surface is always deformed into a prolate spheroid. Drop defor- 
mation toward an oblate spheroid was nevertheless observed in experiments 
Il962h . tively this different t ypes of 



( Allen and Mason , 



shap e, Taylor proposed a theory based on the leaky dielectric model (jTayloii . 
19661 ). which was the b asis for further deve l opments leading to the so-called 



electrohydrodynamics ( Melcher and Taylor . 19691 : Savilk . 1997). A recent in- 
vestig ation of electric deformations on fluid interfaces (jShkadov and Shutov 
2002l ) showed that prolate and oblate shapes in the field direction can be ob- 



tained for both drops and bubbles depending on conductivity and dielectric 
constant ratios. A magnetic field can as well interact with fluid interfaces an d 



goooh . 



form well-organized peak struc tures (IRosensweid . 119971: lLange et al 

or elongate magnetic droplets (IBacri et all 119821 : iBanerjee et al.l . ll999l )~ 

The quantitative description of these deformations also opened new hori- 
zons in fluid metrology by giving the opportunity to characterize, in a non- 
contact way, the mechanical properties of fluid int erfaces, with a particular 
attention to the interfacial tension u sing an electric (Wohlhuter and Ba saranl . 
1992 : Ramos and Castellanos . 19941 ) or magnetic ( Flament et al.l . ll996l ) field. 



A new application domain of these electric and magnetic manifestations is 
digital microflu idics, i.e. the manipulation of two-phase flows and droplets in 
microchannels (Squires and Quake . 20051 ) as for instance, electric actuation 



([Link et alll2006l ) to manipulate flowing droplets. However, in some particu- 



lar cases, the use of electrodes or coils may become diffi cult, particularly for in 
situ characterization in open environments like oceans (ICinbis and Khuri-yakubl . 
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1991) or for single microobjects such as cells (IWottawah et all l200fil ). To 



overcome this difficulty, a second interface deformation approach due to 
radiation pressure effects, has been anticipated. For instance, it is now 
well establish e d tha t acoustic radiation pressure can deform a fluid interface 



(ICinbis et all Il993h and even induce droplet eje ction for drop-on-demand 



devices felrod et all Il989l : iMeacham et~ail 1 20051 ) . Recently, with the emer- 



gence of nanobiotechnologies, this strategy has been successfully extended 
to micrometer and sub-micrometer scales using optical waves. For exam- 
ple, the optical deformation of soft materials in general, and of biological 
micro-objects in particular, brought new insights on t he vi s coelastic proper 



ties and the elast icity of red blood cells (IGuck et all 12000 : iLee et all 12001 



Guck et aL 2001 



The optical radiation pressure of a laser beam focused 
on a fluid interface is also known as an efficient contactless tool for fast 
metrological measurements of inter facial tension ( Mitani and Sakai . 20021 ) 
and viscosity ( Yoshitake et al. . 20051 ) at a microscopic scale, which is partic- 
ularly appealing in the difficult cases of weak interfacial tensions and large 
viscosities. Thus, performing a detailed investigation of the deformation 
of drops by the optical radiation pressure in both the linear and nonlinear 
regimes of deformation, would combine the possibility of extending the fun- 
damental and applied developments of finite volume electrohydrodynamics 
to the optical domain with the advantage of bringing quantitative insights 
on nonlinear deformation of finite size domains by laser waves, a domain 
which is still poorly known. This is the purpose of the present work. 
The optical deformation of liquid inter faces was first observed by Ashkin 
& Dziedzic ( Ashkin and Dziedzic . 19731 ) using a focused laser beam to de- 
form the water/air free surface. They experimentally demonstrated that the 
meniscus always bends toward the medium of lower optical index whatever 
the direction of propagati on of the incident wave. Later on, Zhang & Chang 
( Zhang and Chang . 19881 ) observed the deformation of a water droplet illu- 
minated by a linearly polarized pulsed laser considering two different energies 
(100 mJ and 200 mJ). In these experiments, the optical absorption was as- 
sumed to be small enough to discard thermal effects. At the lower energy, 
Zhang & Chang observed an oscillation of the drop surface due to the prop- 
agation of_capjnary_waves, a phenomenon that wa s further confirmed by La i 
et al. ( Lai et al. . 19891 ) and by Brevik & Kluge ( Brevik and Kluge . 1999). 
In both studies, authors solved the problem numerically using a linear wave 
theory which assumes a linearization of both the flow and optical radiation 
pressure. When the water droplet was illuminated with the highest available 
energy (200 mJ), its front face adopted a sort of conical shape leading even- 
tually to its disruption and the generation of micro-droplets at the tip. In 
order to avoid such a large amount of laser energy and access to stationary 
shapes instead of transient ones due to a pulsed excit ation, new investiga- 
tions were perform e d mo r e recently by Casner et al. (ICasner and Delvilld . 



200 ll : ICasner et aP . 120031 : ICasner and Delvilld . l2003h and Wunenburger et 
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al. ( Wunenbureer et all l2006al |hh using the interface between two liquid 
phases in coexistence close to their critical point. The interfacial tension 
between such phases being extremely small (~ 1CP 7 N.m -1 ) compared to 
the water superficial tension (72 mN.m" 1 ), def ormations of many or d ers of 
magnitude larger t han i n previous experiment ( Ashkin and Dziedzic . 19731 : 
Zhang and Chang . 1988h could be easily observed using a continuous laser 
wave. They experimentally confirmed that an interface always bends toward 
the fluid of lower refractive index whatever the direction of propagation of 
the beam. At large beam powers, typically of the order of 1W , interface de- 
formations become nonlinear. Stable nipple-like shapes were observed when 
the laser wave is incident from the fluid of lowest optical index while propa- 
gation in the opposite direction leads to a needle-like shape, the disruption 
of the interface and a jetting instability driven by the total reflection of 
light within the deformation. As illustrated in Fig. Q], interface disrup- 
tion is still observed in such near-critical two-phase samples during (i) the 
deformation and the subsequent draining of a thin film (of largest refrac- 
tive index) which wets the cell edges near the critical point and (ii) the 
adiabatic growth of a heterogeneous liquid drop (of largest refractive index 
too) during a liquid-liquid phas e transition dr i ven b y a localized composi- 
tion quench induced by a laser foelville et all Il999h . While very different 
in nature, these last manifestations show analogies, at least at the level of 
experimental picture s, with t he de formation and disruption of charged drops 
observed by Zeleny ( Zeleny . 1917 ) and Taylor ( Taylor . 19641 ) under strong 
electrical fields. Indeed, by increasing the electric field, drops deformation 
evolves from rounded to conical shapes emitting eventually a jet o f micro- 
dropl ets at the tip. A theoretical study was proposed by Taylor (ITaylon . 
1964 ) indicating that the stable static solution for a conductive drop sub- 
mitted to both electric and capillary coupling s is near-conical with a semi- 
angle of 49.3°. Since then, many experimen tal ( Odder shedeand Nagell200fi 
Reznik et all. 120041: [Fernandez and Homs"vll2004l:IChen et alll2005h and the- 



oreti c al ( Wohlhuter and Basaranl.ll992l : Ramos and Castellanosl . 1994 ; Saville 
19971 : iFernandez-De-La- Moral . Il992h ~ works have been performed, to analyse 
stability and disruption of these conical shapes as well as on the value of the 
semi-angle v ersus electric prop erties of the liquids (see for instance Fig. 2 in 



Stone et al. (IStone et al.l.ll998f) for the dielectric c ase and the review of Fer 



nandez de la Mora ( Fernandez-De-La-Mora . 19921 ) for conducting droplets) 
Note finally, that deformed interfaces with conical shape can be found in the 
absence of electromagnetic excitation. A first example is the near-conical 
shapes observed during the deformation of the interface between immiscible 
fluids by selective withdra wal, with a straw who s e tip is suspended abov e 



20071 : ICohen and Na.gell. [2002) 



the unperturbed interface (ICase and Nage 
Other fluid dynamics examples are drop brea kup from a nozzle (lEggersl. 1997) 
or sin k flows in the presence of an interface ( S.CourrechduPont and Eggersl . 



2006), thus illustrating the emergence of a sort of robust and general topolog 
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ical transition ( Zhang . 20041 ) of fluid interface deformations under localized 



forcing. Motivated by some amazing effects of the radiation pressure illus- 
trated in Fig. Q] and suspecting that finite size effects may promote the 
observed near-conical shapes, we propose a numerical investigation of op- 
tical deformation of sessile liquid drops by continuous laser waves. Indeed, 
beyond the simple deformation of spherical interfaces, we demonstrate in the 
present work that optical stretching can lead to near-conical shapes. Above 
a radiation pressure threshold, drop disruption is also observed numerically, 
in the absence of total reflection of light however, demonstrating a new phe- 
nomenon of interface instability. The optical analogue of oblate deformation 
by an electric field, induced here by an optical squeezing, is investigated as 
well. At large radiation pressures, the local squeezing leads to the forma- 
tion of torus-like shapes. Section II is devoted to the physical model used 
to predict the deformation of drops by the optical radiation pressure. Sec- 
tion III briefly summarizes the numerical algorithm used here, based on the 
Boundary Integral Element Method (BIEM). Results on drop deformations 
for both stretching and squeezing are presented and discussed in Section IV. 

2 Physical model 

The droplet configuration under consideration, together with the notations 
used throughout this work, are represented in Fig. [2j 

To describe the hydrodynamics of the laser/fluid interaction, let us first 
consider a Gaussian continuous laser wave, of beam waist uq, and make di- 
mensionless all lengths involved in the problem using uq. This laser wave 
is supposed to impinge on the interface between a semi-spherical dielectric 
viscous drop (marked as fluid 1) of dimensionless radius a (a = where Rd 
is the drop radius, a = 2 in Fig. [2]) and a surrounding dielectric viscous fluid 
(marked as fluid 2) of dimensionless horizontal and vertical extensions R and 
H. We assume that the three-phase contact line is motionless. Since the in- 
tensity distribution of the laser beam, centered on the drop, is supposed to 
be axisymmetric as in usual situations, cylindrical coordinates (e r , e z , e<p) 
with their origin O located at the center of the drop are used throughout this 
work so that any point x is located by (r, z, 4>) in this reference frame of co- 
ordinates. Governing equations are written in a dimensionless form using (i) 
the laser waist too as the characteristic length scale, (ii) the viscous relaxation 
velocity u* = -fa as the reference velocity associated to the characteristic 

timescale t* = < ^ UJ ° and, (iii) the reference pressure p* = , i = 1,2 
to rescale the pressure pi in each phase. Here, 7 is the interfacial tension 
between fluids 1 and 2, )Uj is the dynamic viscosity of fluid i and (y) = Ml "^ M2 
is the average viscosity of the fluids. It could be pointed out that, depending 
on the fluids and/or the large beam powers under consideration, (those given 
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(a) 

N 2 <N 1 



111 II I 




Figure 1: (a) Dynamics of the draining and the deformation of a thin wet- 
ting film by the optical radiation pressure of a continuous Ar + laser beam. 
The initially flat interface (not shown) takes a rounded shape which further 
becomes near-conical and eventually destabilizes to give birth to a jet. The 
fluids are phase-separated l iquids close to a liquid-liquid critical point (see 
( Wunenbureer et all l2006aT l 'or a detailed description of the system). The 



thin film results from the wetting transition occurring close to the critical 
point. Time delays from the first image are respectively 12 s, 18 s, 22 s, 23 s 
and 24 s (from top left to down right). Control parameters are P = 590 mW, 
ujq = 4.2 /im and T — Tq = 5 K, where Tq is the critical temperature. The 
arrow indicates the direction of propagation of the exciting beam, not seen by 
using a colored filter, (b) Dynamics of a heterogeneous drop growing adiabat- 
ically after an optical quench in composition inside the coexi stence curve of 



a bin ary liquid mixture driven by a continuous Ar + laser (see (|Delville et al 



19991 ) for experimental details). The drop is heterogeneously nucleated on 



the bottom cell window and grows inside the beam. During the late-stage 
growth, its shape shifts gently from hemispherical to near-conical until desta- 
bilization of the interface by the radiation pressure. Time delays from the 
first image are respectively 202s, 266s and 312s (from left to right). Control 
parameters are P = 880 mW, u>o = 12-4 /im and T = 293.6 K. The observed 
laser beam propagates upwards. 



in Fig. 1, for instance), temperature effects may disturb or even overcome 
the mechanical effect of light on fluid interfaces, thus making a coupled heat 
and momentum transfer description necessary. The first expected additional 
effect is a direct laser heating due to the optical absorption of the drop 
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Figure 2: (Color online) Schematic representation of a captive drop (fluid 1) 
immersed in a second liquid (fluid 2) and submitted to the optical radiation 
pressure of a laser beam centered on the drop axis and initially focused at 
the interface. Sci, Sc2, and Sj respectively denote the solid boundary with 
fluid 1, fluid 2 and the interface between the two fluids. The initial drop is 
hemispheric. The three-phase contact line is assumeed to be pinned. Lengths 
are made dimensionless using the beam waist u>q. 



and/or of the surrounding fluid. Fluids must indeed be transparent at the 
used optical wavelength as it is generally the case for classical liquids in the 
visible window. Typically, the optical absorpt ion of water, used in Ashkin 
Sz Dziedzic and Zhang h Chang experiments Ashkin and Dziedzic ( 19731 ): 



Zhang and Chang ( 19881 ). is of the order of 3.10 3 cm 1 in the visible region, 



while that of the micellar phases used in the examples illustrated in Fig. 1 is 
3.10 _4 cm _1 . In the latter, the ov erheating induced by a beam power of the 



order of 1 W is smaller than 0.1 K (IChraihi et all 12007ft . Direct laser heating 



effects can then be discarded, even in critical fluids as far as temperature is 
not too close to the critical one. Thus, we consider in the following all liquid 
properties (7, pi, pi, N; L ) as constant in the presence of laser light, pi and 
Ni being respectively the density and optical refraction index of fluid i. A 
second coupling is the thermocapillary effect. Since the interfacial tension 
7 is a function of the temperature, local laser heating may drive interfacial 
tension gradients inducing stresses on the drop interface and its subsequent 

for classi- 



deformation (lLoulergue et al.l . ll98ll ). A typical value of I7 



K dT J 



cal fluids is 10 3 K 1 Sammarco and Burns ( 19991 ). which leads to negligible 



thermocapillary effects considering the above mentioned laser overheating. 
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For the case illustrated in Fig. 1, we already found th at thermocapil l ary in - 
terface deformation is negligible near the critical point Chrai'bi et al. ( 20071 ). 
Consequently, we can safely discard temperature effects without affecting 
the generality of our purpose. Finally, when getting very close to the critical 



point, the capillary length Iq 



where Ap is the density contrast 



and g the acceleration due to gravity, vanishes while thermal fluctuations 

Aarts et al.l i2Q()4 ). 



increase. Then, the interface roughness It 



lk B T 



1 



where /cgT is the thermal energy, increases too an d may dominat e capil lary 
effects. Using the critical data given in Reference IChrai'bi et al.1 (120071 ) for 
deformations presented in Fig. 1, Iq >> It requires T — Tq » 3.10 _2 -f£T, 
a condition which has always been fulfilled in experiments. Conversely, for 
classical fluids, 7 lies between 10 and WOmN/m and Ap is about 10 2 to 
10 3 kg /m 3 . As a consequence, It is orders of magnitude smaller than Iq so 
that thermal fluctuations are also negligible for this class of fluids. 

In addition, we assume in this study that inertial and gravity effects 
are negligible at the micrometric scale, which implies that the Reynolds 
and Bond number are small compared to unity. Along with the condition 
lc » It which is also automatically satisfied as indicated above, we consider 
an incompressible quasi-static Stokes flow in each phase. 

Therefore, the hydrodynamics of each liquid phase is described by the 
Stokes and mass conservation equations, respectively given by: 







and 



- Vft + V 2 Ui 



V • u ; = 



1,2, 



1,2, 



(1) 



(2) 



where u, and qi are respectively the dimensionless velocity and the pseudo- 
pressure in fluid i. This pseu do-pressure contains t he ele ctrostrictive contri- 
bution of light in dielectrics ( Landau and Lifshitz . 1960l ) which was demon- 
strated to have no incidence on the shape or the height of the deform ed in- 
terface iLai et al.1 . flQSfll : iBrevik and Klugd . ll999l : IChraibi et allEuOTl ). This 
pseudo-pressure is defined as 



P* 1 



€ °F 2 n 



(3) 



where ej = Nf represents the relative dielectric permittivity of fluid i, eo is 
the permittivity of vacuum and Ef is the quadratic magnitude of the electric 
field in fluid i averaged over an optical period. Considering the classical 
expression of the divergence free hydrodynamic stress tensor Tj, 



T ; = - qi I + (Vui +* Vui) 



(4) 
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we can write the boundary condition at the interface as follows: 
2 



-(AT] 



n 



n • n 



n(r) — II(r), 



(5) 



1 + A 

where n is the unit vector normal to the interface directed from fluid 1 to 
fluid 2, and A = ^ is the viscosity ratio. Eq. ([5|) simply expresses the 
fact that the normal stress on the interface is balanced by capillary forces 
and the optical radiation pressure, respectively represented in dimensionless 
forms by k(t) and II(r). Indeed, in the right hand side of Eq. ([5]), k{t) 
represents the dimensionless double mean curvature of the interface, given 
by: 



Id 

r ar 



rz 



(6) 



where z' = 4p is the local slope of the interface. The second term, n(r), is 
the contribution of the dimensionless optical radiation pressure at the drop 
interface. As the interface deformation is experimentally found to be directed 
toward the fluid of smallest r e fractive index whatever the direc tion of propa- 



gation (lAshkin and Dziedzid . Il973l : ICasner and Delvilld . l200ll ) , we consider 



the Minkowski point of view which states that the photon momentum in a 



diele ctric medium varies linearly with the refractive index (| Campbell et al 



for a review on the Abraham -Minkowski controversy see for instance 
the review of Brevik (iBrevikl Il97flh . The amplitude of U(r) is neverthe- 
less affected by the direction of propagation through its dependence on the 
incidence and transmitted angles. In addition, II(r) also depends on the 
polarization of the laser wave. In the present work, we preserve the axial 
symmetry of the laser/fluid interaction by assuming a circular polarization 
of the laser wave. Finally, for the sake of simplicity, we assume refringence 
of light at the drop interface at any incidence angle. By doing this, we elim- 
inate situations where total reflection of light may occur at the interface as 
those presented in Fig. Q]for example. Therefore, the laser wave is supposed 
to propagate from the optically less dense fluid. Note that this choice is not 
restrictive at all; by removing this assumption, we would find the same type 
of drop deformations but with an asymmetry in amplitude due to the non- 
linear behavior of the transmission and re fl ection Fresnel coeffi cients with 
the angle of incidence ( Casner et al. . 2003 : Chrai'bi et al. . 200?! ). Another 



reason for this choice is made clear in Sec. IV. 2 where a new type of laser- 
induced interface instability is advanced, the one demonstrated previously 
being triggered by total reflection of light. 

On the one hand, to stretch a drop by light, and thus deform the interface 
outward (Sec. IV. 1), the wave propagates downward. In this case, one has 
-Ai > N2, and the dimensionless expression of the optical radiation pressure 



is nw = n 



U dn {r) 



dn 



(Of 



where IL (r) is given by: 



I(r) 



cos 6i(2N 2 cos 6i - ip dn {N 2 cos 0* + Ai cos 9 t )). (7) 
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On the other hand, to squeeze a drop and thus deform the interface inward 
(Sec. IV. 3) we choose N\ < N 2 and consider a beam propagating upward. 
In this case the expression of the rescaled optical radiation pressure becomes 
II(r) = U u P(r)^ where U u P(r) is given by: 



n" p (r) = -W- cos 0i(2JVi cos Q, L - ^ up (iVi cos t + N 2 cos OA). (8) 
c 

In Eqs ([THE]), we have respectively denoted by Q% and 6 t the incidence and 
transmission angles; c is the light celerity in vacuum. One has Oi = arctan(V) 
while t depends on the direction of propagation. One has Ot = arcsin(r/ sin Oi) 
or t = arcsin(isin^j) respectively for upward and downward propagations, 
where 77 = N\/Ni is the refractive index ratio. Moreover, I(r) represents the 
intensity of the Gaussian laser beam. By neglecting the weak z-dependence 
of the beam radius, its expression is given by: 

op , 

J(r,*)«J(r) = — e~ 2r , (9) 

where P is the beam power. In addition, ty up and ty down are the Presnel 
transmission coefficients of energy fluxes. These coefficients are obtained 
from the ratio of the transmitted to incident normal components of Poynt- 
ing's vectors and are expressed as 

yuup = yudown = 2JVijV 2 COS 0j COS t 2jVijV 2 COS Oj COS t 

(iVi cos Oi + N 2 cos t ) 2 (N 2 cos { + N x cos t ) 2 1 ' 

for a circulary polarized beam. In order to quantify effects of the laser wave 
on the interface deformation, it is convenient to define the electromagnetic 
to Laplace pressure ratio. This ratio £, taken at r = {Oi = Ot = 0), is 
defined as: 

C= n(r = 0) = , t = 1,2, 11 

7TCW07 (jV2 + JVi) 

where Ni refers to the optical index of the incidence fluid. When no slip is 
assumed at the interface, along with the fact that fluids are immiscible, it 
follows that the velocity u, of the interface Sj is equal to that of each fluid 
particle on Si, i.e.: 

u(x) = m(x) = u 2 (x) xeS/ (12) 

Moreover, the movement of the interface is described using a Lagrangian 
approach. It consists in following each fluid particle of the interface in its 
Lagrangian motion according to the kinematic condition: 

r/x 

— = u(x) xGSj. (13) 
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This condition indicates that the interface is advected along with the flow 
until the equilibrium is reached for which normal velocities along the interface 
are zero, i.e. u(x).n = 0, x € Si. Finally, we assume a classical no-slip 
boundary condition on all the solid boundaries of the domain 

Ui (x)=0 x = 1,2. (14) 

The above system of equations is solved using a Boundary Integral Element 
Method (BIEM). Due to the axial symmetry of the laser/drop interaction, 
it consists in an axisymmetric integral formulation making use of the fun- 
damental solution of Stokes equations. The solution is sought with a con- 
stant boundary elements discretization technique according to the numerical 
scheme described below. 



3 Numerical algorithm 



A brief description of the numerical algorithm is presented in this section. 
For more extensive details on the BIEM applied to a two-phase axisymmetric 
flow, the reader may refer to the review by Tanzosh et al. on the solutio n 



of free surface flow problems using this technique ( Tanzosh et al. . 19921 ) 



The BIEM reveals to be an excellent tool to solve interfacial flow problems 
with high resolutio n as reported in the analysis of flow involving electric 
and magnetic fields (Sherwoodl . 1987 ) or buoyancy ( Manga and Stone . 19941 ) 
JKoch and Kochl . Il994h . 

Because solutions to Stokes equation can be formulated in terms of Green's 
functions, we can rewrite the governing equations as a system of integral 
equations over the boundaries of the computational domain. When doing 
so, the boundary integral form of the Stok es equation for fluid i (i = 1,2) 
can be written as follows (jPozrikidisl . Il992h 



1 

-u, 

2 



x 



U • (Ti-n^dSy - / rii-K • mdS y , (15) 
'Si+Sa J Si 

rij being the unit normal vector directed toward the outside of fluid domain i. 
In this last expression, U and K are second and third order tensors forming 
the Green's kernel for velocity and stress ass ociated to the St okes equation. 
These two tensors are respectively given by (IPozrikidisl . ll992h : 



U(d) 
K(d) 



1 .1 dcL 
3 ddd 



(16) 
(17) 



In these two relationships, d = x — y, y(r y , z y ) is the integration point. Once 
boundary conditions on Si, Sci and Sc2 are used, the two-phase Stokes 
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problem can be written in the following compact form: 

u(x)= I \J.n(K(r y )-n(r y ))dS y + 
J Si 

— ^— [(1-A) / n.K.udS y + x[ U-(Ti-n)d5j, 

1 + A J Si JSci 

I V ■ (T 2 ■ n)dS y ]. 

J Sc2 



[11 



Here, 5/, Sex and 5c2 must be understood as axisymmetric surfaces instead 
of their trace in the plane of Fig. [2J In Eq. (|18|) . the first term on the right 
hand side describes the flow contribution from surface tension and radiation 
pressure, whereas the second term accounts for the shear rate contrast at the 
interface. As expected, this second term vanishes when there is no viscosity 
contrast between the two phases (A = 1). The third and fourth terms account 
for the shear occurring on 5c l and 5c2 as a result of the no-slip boundary 
condition. 

The solution of the problem consists in the computation of the velocity, u, on 
the interface as well as the stress over all the boundaries Sj,Sci and 5c2- This 
is performed once Eq. (JTHJ) has been discretized using boundary elements. 
Here, we use constant boundary elements for which pressure, velocity, and 
hence stress, take constant values on each element, equal to that at the 
central node. The overall numerical procedure can be summarized as follows. 
(1) All boundaries Sj, 5ci and 5c2 are discretized with line segments of total 
number N. The fluid-fluid interface 5/ is parameterized in terms of an arc 
length, s, so that the double mean curvature can be accurately computed 
using the following expression: 



dt 

ds 



+ 



z'(r) 



r(l + z(r)' 2 ) 1/2 ' 



(19) 



where t is the tangential vector to 5/. The number of mesh points on 5/ 
is 70 for a typical computation. Each horizontal and vertical solid bound- 
aries are meshed using approximately 50 uniformly distributed points. An 
increase in the mesh resolution for the interface and solid boundaries did 
not show any significant change in the results. (2) At each time step, the 
azimuthal integration of the integrals arising once Eq. (18) is discretized 
with N constant bound ary elements is performed analytically ( Lee and Leall . 
19821 ). ( Graziani . 19891 ) reducing Eq. (Tl8l) to line integrations which a re fi- 
nally performed using Gauss quadratures (iDavis and Rabinowitzl . 1 1 9 H). El- 
liptic integrals resulting from the azi muthal integration are evaluated using 
power series expansions (lBa,krl . ll985h . (3) Once all integrals are calculated, 
the linear system Aw = b is formed where w is the vector of the unknown 
boundary values of velocities on Si and normal stress and shear on 5/, 5ci 
and 5c2- The matrix A results from boundary integrals computed with the 
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procedure detailed above. The vector b is obtained from the product of 
boundary integrals and the known part of the normal stress jump on Sj as 
it appears in the first term in the right hand side of Eq. (jTHJ) . This linear 
system is solved at each time step using a direct method based on a LU 
decomposition. This step provides u(x) on Sj. (4) Finally, the motion of the 
interface is captured using the kinematic condition (Eq. fl~3j) and an explicit 
first-order Euler time scheme the discrete form of which can be written as 
follows: 

r{t + At) = r(t) + u r (t)At, (20) 
z(t + At) = z{t) + u z (t)At, (21) 

where At is the time step. This time step is typically chosen to be about 
20 times smaller than t*. Once the interface has been moved (stretched or 
squeezed), the mesh is reconstructed using a smoothing procedure with local 
cubic splines. Starting from an initial semicircular interface at rest, the laser 
beam is switched on at t = and the interface starts to deform towards 
the fluid of smallest refractive index. The overall algorithm is repeated until 
final equilibrium is reached and this is achieved when 

u(x).n < e x € Si, (22) 

where e << 1 is a user fixed parameter. 

4 Results 

In this section, we first present results on drop stretching by the optical 
radiation pressure. An outward deformation of the liquid drop is obtained 
when its optical index of refraction is larger than that of the surrounding fluid 
(i.e. Ni > N2), resulting in the formation of near-conical shapes at different 
dimensionless initial drop radii a and refractive index ratios, characteristic 
of most dielectric fluids. An instability arising from the destabilization of 
the drop above a threshold value of £ is also demonstrated. Results on drop 
squeezing corresponding to Aq < A2 are presented in a second part. In this 
case, drop deformation becomes concave and may reach a torus-like shape. 

4.1 Drop stretching by the optical radiation pressure 

Simulations of outward deformations were performed with A = ^ = 2 and 
rj = jf* = 1.1. This value of the index ratio, corresponding for example to 
a toluene droplet immersed in water, was chosen as a representative couple 
for the water/oil interface. The generality of the investigation is nonetheless 
totally preserved because A has no influence on the stationary shape of de- 
formed drop (see the left Inset of Fig. [3]) and no bifurcation has been observed 
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in our study. Variations in 77 just lead to qualitative differences in optical 
radiation pressure effects, as illustrated in the near-conical shape study (see 
Sec. IV. 2). Stationary shapes of an optically stretched liquid drop of dimen- 
sionless radius a = 0.5 are illustrated in Fig. [3]for different optical to Laplace 
pressure ratios £. In Fig. 02 one clearly see that the shape of the drop pro- 
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Figure 3: (Color online) Optical stretching of a drop of dimensionless radius 
a = ^ = 0.5 (R(i is the drop radius, uiq the beam waist). The number 
associated to each outward deformation refers to the pressure ratio £. The 
initial spherical shape is represented by the dashed line. The arrow indicates 
the direction of propagation of the laser beam. Left Inset: steady shapes 
at £ = 50 for two viscosity ratios showing no effect of A on the steady-state 
solution. Right Inset: variation of the optical radiation pressure rescaled by 
its value at normal incidence, versus the incidence angle for ^ = 1.1. 

gressively varies from rounded (£=10, 20, 40) to near-conical (£=500). These 
stationary shapes result from the competition between capillary forces and 
optical radiation effects. Increasing £ increases the height of the deformation 
and thus induces in turn an increase of the curvature of the interface. This 
effect is maximized at the tip of the deformation due to the Gaussian profile 
of the light intensity and the normal incidence at the interface at r = 0, as 
illustrated in the right inset of Fig. [3] . Therefore, an important increase in 
£ eventually leads to pointed drop shapes, assuming the fixed contact line 
hypothesis and the finiteness of the drop volume. Although the stress on the 
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interface is here of optical origin rather than electrically or magnetically in 
duced, the deformation of drops by the optical radiation pressure shows strik- 
i ng similarities with thos e obse r ved in electro- and magneto-hydrodynamics 
(Wo hlhuter and Basaranl . 19921 : Ramos and Castellanos . 1994 : Stone et al 



19981 : ISherwoodl . ll987h . There are nevertheless two major differences. First 



the relative extension of the exciting field versus the drop size, considered 
as infinite in the case of electric or magnetic fields, is intrinsically finite due 
to the Gaussian shape of the laser intensity in the current analysis. More- 
over, the associated radiation pressure nonlinearly decreases with the inci- 
dent angle through the Fresnel transmission coefficient, drastically reducing 
the mechanical effect of the beam at the edge of the deformation when its 
aspect ratio increases significantly (see the right Inset of Fig. \3§. Second, the 
pointed shape observed on the axis does not originate from a local increase 
of the electric field at the tip as in the electrostatic case because at a given 
£, the optical radiation pressure at the tip, where the incidence angle is close 
to zero, remains regular. 
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Figure 4: (Color online) Optical stretching of a drop of dimensionless radius 
a = 1 (left) and a = 2 (right). The number associated to each outward 
deformation refers to the pressure ratio £. The initial spherical shape is rep- 
resented by the dashed line. The arrow indicates the direction of propagation 
of the laser beam. 



Fig. H] shows the outward deformation of a drop with dimensionless 
radius a = 1 and a = 2. Qualitatively, droplet deformations are very similar 
to those of Fig. [3l Increasing the dimensionless radius of the drop leads to a 
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decrease of the £ value required to achieve a given dimensionless deformation 
height. In other words, a larger drops will have a higher aspect ratio h/a 
than a smaller one for the same £, h being the equilibrium height defined 
as h = z(r = 0, t — ► oo). The £ variation of the drop aspect ratio h/a, is 
illustrated in Fig. [5]for different drop radii a. The fact that the h/a(£) curves 
depend on a can be explained by the distribution of radiation pressure on 
the drop, as shown in the inset of Fig. When oCl, the electromagnetic 
intensity applied to the drop is almost uniform. The radiation pressure 
mainly depends on the incidence angle which vary for an initial spherical 
drop from 0° at the tip to 90° at the contact line. So, it is nearly uniform 
except near the contact line where it decreases down to zero, as shown in 
the inset of Fig. El The associated electromagnetic normal stress is almost 
uniform all along the drop surface. Conversely, when a>l, the radiation 
pressure mainly depends on the gaussian intensity of the laser beam, as 
the incidence angles near the tip remain close to zero. The electromagnetic 
normal stress is localized on the top of the drop. Consequently, the drop is 
more deformed for a 3> 1 than for a <C 1. 




Pressure ratio £, 



Figure 5: (Color online) Variation of the reduced drop stretching h/a versus 
the pressure ratio £ for different dimensionless drop radii a. The inset shows 
the variation of the optical radiation pressure rescaled by its value at normal 
incidence versus r/a. 
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4.2 Optical cone formation and interface instability 

As illustrated in Figs. [3] and HI drop deformations with near-conical shape 
seem to emerge when increasing the pressure ratio £ whatever the beam 
radius. It could be objected that such a generic shape requires an optical 
radiation pressure behaving as 1/r to balance the Laplace pressure, while 
the optical coupling shows a dependence in both exp(— r 2 ) and a nonlin- 
ear function of the incident angle. However, as experimentally illustrated 
in Fig. Q] and in the following, under certain conditions, conical shapes 
can approximate the true deformation profile with a high degree of accu- 
racy. The existence of such type of deformation is still surprising if we 
take into account predict i ons for dielectric drops subjec t ed to electric fields 
( Wohlhuter and Basaranl . 19921 : Ramos and Castellanos . 1994 : Stone et al 



Sherwood! 19871 ) which indicate that conical shapes cannot exist for 



dielectric constant ratio smaller than 20. Here, the situation is nevertheless 
different since the beam profile is inhomogeneous and the deformation ef- 
ficiency depends on the beam extension. To quantitatively investigate the 
emergence and the existence of cone-like shapes, we implemented the follow- 
ing procedure. For a cone, the semi-angle is defined as 90 — arctan(|z'|) (in °), 
where z' is the slope of the generating line. Considering a drop of dimension- 
less radius a, we plot 90 — arctan(|z'|) for increasing £. A cone-tip signature 
appears through the emergence of a plateau, corresponding to a constant 
local slope z' along the interface, which defines its semi-angle. As illustrated 
in Fig. [6] for a = 0.5 and a = 2, i.e. for beam waist respectively larger and 
smaller than the droplet radius, a plateau emerges at a threshold £p ~ 50 
(a = 0.5) and £p ~ 8 (a = 2) for N1/N2 = 1.1. Its wideness increases 
for £ > £p while the corresponding semi-angle remains almost constant. 
This saturation of the semi-angle at large £, observed at any investigated 
dimensionless drop radius, demonstrates the robustness of the near-conical 
deformation. Results in Fig. [6] also suggest that the value of the semi-angle 
varies with the dimensionless drop size a. In Fig. we present stationary 
deformations of liquid drops of different dimensionless radii when £ > £p. A 
near-conical shape is systematically observed. These results also illustrate, 
in a different representation, the variation of the vertical semi-angle of the 
cone versus the radius of the drop. To be more quantitative, Fig. [8] shows 
that the semi-angle (i) remains finite at large values of £ whatever the value 
of a and (ii) exhibits two different asymptotic values at large and at small a, 
the larger one, at small a, logically corresponding to the steep variation of 
the mean slope of the deformation h/a{^) shown in Fig. El and the smaller 
one at large a to the smooth variation of h/a(£). For small a, the asymptotic 
value is roughly 26° while it is close to 5° at large values of a. 

The finiteness of the asymptotic value of the cone semi-angle at large £ 
is the most striking behavior of this class of laser induced interface deforma- 
tions. The fact that the optical radiation pressure pulls the cone tip with a 
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diverging strength when £ — ► oo would make one indeed think that the cone 
slope also diverges. This behavior has a consequence on the cone stability. In 
fact, above a second threshold value £/ of £, the drop deformation becomes 
unstable leading to an interface breakup, as illustrated and discussed below. 

In the a <C 1 limit, the experimental conditions used to deform the 
drop are similar to those used with electric fields in capacitors since the 
radial extension of the exciting laser wave is significantly larger than the 
drop size. Moreover, the asymptotic value of the semi-angle 26° is close to 
the minimum one, 30°, found for dielectrics in electric fields. This result is 
surprising since the dielectric constant ratio required to reach this minimum 
angle is (N1/N2) 2 = 1.21 while in the electrostatic case it corresponds to 
17.6. In addition, the asymptotic value of the semi-angle in the case of optical 
excitation at small a was found to remain equal to 26° for {N 1 /N 2 ) 2 = 17.6, 
an unrealistic optical situation for dielectric liquids. This mismatch between 
electrical and optical excitation is likely due to the fact that the mechanisms 
involved in the formation of near-conical shapes are different. However, even 
if we do not have yet any quantitative explanation, these observations, at 
least at the level of experimental observations, suggest the existence of a 
possible common phenomenology of the deformation of dielectric drops by 
electric and optical excitation. To give another insight on this appealing 
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Figure 7: (Color online) Stationary tip shapes at N1/N2 = 1.1 for different 
dimensionless radii at £ values above the semi-angle threshold emergence £p. 



comparison, we analyzed the index contrast variation of the semi-angle for 
different dimensionless drop radii (Inset of Fig. [8]). The refractive index ratio 
was varied from 1.05 to 2, this last alue corresponding to an upper limit for 
the free surface of dielectric liquids (the largest index of refraction we know is 
that of diiodomethane CH2I2, N = 1.7425 at 20°C). Given a dimensionless 
radius a, the measured angle weakly increases with N\/Ni when a = 1, while 
it is almost constant for a « 1. This last behavior can be explained by the 
two following facts: (i) the laser incidence experiences all angles between 0° 
and 90° when the beam waist is larger than the drop radius and (ii) the light 
intensity is almost uniform over the drop when a « 1. 

All the near-conical shapes presented above were obtained before the 
previously mentioned instability threshold. They were thus stable. Given 
a dimensionless drop radius a, when the pressure ratio £ reaches what we 
called the second threshold value £/ (the first one, £p, was attached to the 
emergence of a near-conical shape), a destabilization of the deformation is 
observed, leading to a disruption of the interface at the tip, as illustrated 
in Fig. [8l A refinement in mesh and time steps did not show any change 
on this behavior. Even if the optical radiation pressure is known to be able 



to destabilize a flu i d inter face and form a liquid jet (ICasner et all 2003 



Wunenbureer et all l2006bh . this numerical result is quite different because 
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Figure 8: (Color online) Semi-logarithmic variation of the semi-angle of near- 
conical stationary shapes versus the dimensionless drop radius a (rj = 1.1). 
Inset: variations of the semi-angle versus the refractive index ratio over the 
range of investigated dimensionless drop radii. 



jetting was previously induced by total reflection of light within the defor- 
mation when the laser beam propagates from the fluid having the largest 
refractive index towards that of lowest refractive index, as it is the case 
in Fig. [H To prevent this effect, here the opposite situation was chosen. 
The beam propagates from the fluid having the lowest refractive index to- 
wards that of largest refractive index, a situation where "nipple-like" inter- 
face defor mations are generated at la rge beam powers for fluids of infinite 
extensions ( Casner and Delville . 20031 ) and for which further increase of the 
power keeps the nipple like shape stable. The main difference with the 
present investigation is that now, the volume of one of the two fluids, the 
drop, becomes significantly small compared to the other one. Increasing £ 
still increases the height of the deformation, but the drop turns to adopt a 
near-conical steady shape, the semi-angle being surprisingly preserved along 
a wide plateau which extension increases towards smaller values of r while 
increasing £ (see Fig. [6]). This complex balance between the Laplace pres- 
sure and the laser incidence dependent radiation pressure is preserved up 
to £j, as illustrated in Fig. [9) The saturation of the semi-angle at large £ 
observed numerically implies that the radiation pressure applying along the 
cone slope linearly increases with £ whereas the cone shape does not evolve 
anymore: this qualitatively explains why a destabilization of the strongly 
deformed drop occurs beyond the threshold value £/, as illustrated in Fig. [9l 
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The quantitative investigation of this drop disruption under a very high op- 
tical radiation pressure regime and the subsequent disruption, which clearly 
deserve a devoted study, will be presented in future development. 
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Figure 9: Dynamics of the disruption of the near-conical shape of a drop of 
dimensionless radius a = 0.5 at different times (the reduced time step is 0.02). 
The lowest profile is the stable stationary shape obtained for £ = 500. At 
£ = 700, the destabilization of the deformed drop occurs. The corresponding 
profiles (from down to top) have been translated for clarity. The arrow 
indicates the direction of propagation of the laser beam. 



4.3 Drop squeezing by the optical radiation pressure 

In this section, we investigate the deformation of a liquid drop of refractive 
index lower than that of the surrounding liquid. We preserved A = ^ = 2 for 

calculations and reverse the refractive index contrast by taking rj = ^ = 0.9 
in order to force a drop squeezing instead of stretching. Moreover, we re- 
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stricted our investigation to the case a > 1 in order to avoid any significant 
effect of multiple reflection inside the droplet. As before, the generality of 
the purpose is preserved because the value of A has no influence on the 
stationary shape of deformed drops and variations in r\ just modify quanti- 
tatively the optical radiation pressure effects. Fig. [TO] shows the resulting 
variations of the stationary shape of a liquid drop of dimensionless radius 
a = 1 and a = 3 at various pressure ratio £. The spherical shape flattens at 
low values of £. For larger beam power, the curvature reverts and the shape 
becomes concave. At large £, the concavity reaches the solid boundary and 
flattens more and more in its central region. As in the stretching case, the 




Figure 10: (Color online) Optical squeezing of drops of dimensionless radius 
a = 1 and a = 3. The number associated to each inward deformation refers 
to the pressure ratio £. The initial spherical shape is represented with a 
dashed line. The arrow indicates the direction of propagation of the laser 
beam. 



main difference between the deformations obtained at different dimensionless 
drop radii is that the amplitude of the deformation at a given £ decreases 
with the drop radius and reaches some asymptotic behavior. In Fig. [TT1 we 
have represented the evolution of the reduced height of the drop h/a as a 
function of the pressure ratio £ for drop radii 1 < a < 3. For a > 2, the 
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Figure 11: (Color online) Variation of the reduced drop squeezing h/a versus 
the pressure ratio ^ for different dimensionless drop radii a. 



evolution of the reduced height shows little change while increasing a. A 
behavior similar to that was already observed for drop stretching. These 
observations can as well be explained by the decrease of the incidence angle 
range of the illuminated area of the drop while increasing its radius. For 
a = 1, larger radiation pressure is required to obtain the same amplitude of 
the interface deformation. The former argument proposed for drop stretch- 
ing can be used again to explain the observed increase in £ required to reach 
equilibrium: when a >> 1, the local illumination of the drop leads to larger 
deformations compared to the case where the field is more uniform (a = 1). 
While optical stretching leads to near-conical shapes at large £, the squeez- 
ing of a drop should give birth to the formation of a stable "optical torus". 
Indeed, by assuming a fixed contact line on the substrate, we implicitly pre- 
vent drop spreading, thus promoting the formation of an annular rim. The 
film thinning at the center should as well induce local "de wetting" when th e 



deformation reaches the substrate, as in the electric case (jYeo et all 120071 ). 
Here, we assume that molecular forces are much larger than any external 
forcing at a nanometric scale and, speculate that Van der Waals f orces will 
dominate the final stage of the deformation ( Seemann et al. . 200ll ). a result 
which could be retrieved by including the disjoinin g pressure i nto th e Stokes 
equations or boundary conditions at the interface ( Yeo et al. . 200?! ). Conse- 
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quently, even if the formation of stable torus from radiation pressure effects 
is still speculative and deserves an experimental demonstration, the present 
investigation shows that lasers allow for drop squeezing with morphologies 
which were, to the best of our knowledge, never even suggested while using 
electric fields. 



5 Conclusions 



Although many appealing applications of optically induced drop deforma- 
tions can be advanced, a mong them, the op tical stretcher tool which was 

to deform red bl ood c ells and 
2004h or the 



developed by Guck et al. (IGuck et all 120011 



discriminate between sane and cancerous cells (Lincoln et al 



contactless viscoelastic micro-characterization of fluids dSakai et al.. 2001 



very few theoretical or numerical studies were performed in this field, espe- 
cially in the nonlinear regime of deformation which is strongly influenced by 
finite-volume effects. The objective of this paper was thus to provide new 
elements for understanding such type of drop deformation and go even fur- 
ther in order to illustrate the specificities of drop deformation by the optical 
radiation pressure. Both deformation cases, the prolate one when the drop 
is stretched and the oblate one when it is squeezed, were studied at steady 
state as a function of the amplitude of the optical radiation pressure, normal- 
ized by the Laplace pressure, and for varying drop radii. We found that the 
elongation of stretched drops significantly varies with the beam waist and 
can adopt a near-conical shape at large optical radiation pressures, as sug- 
gested by the experimental illustrations of Fig. [TJ Contrary to the classical 
electrodynamics case, where a minimal dielectric constant ratio is required 



to reach cone shapes ( Stone et al. . 19981 ). these shapes are observable at any 



optical refractive index ratio. The semi-angle was found to be a decreas- 
ing function of the drop radius, showing two asymptotic values. When the 
drop is much smaller than the beam, the semi-angle is close to that obtained 
for electrically deformed dielectric drops. Above a threshold in radiation to 
Laplace pressure ratio, a disruption of the drop is observed. This behavior is 
appealing because even if a si milar phenomenon was alr eady observed exper- 
ime ntally by Zhang fc C hang ( Zhang and Chang . 1988h on drops, Casner et 



al. (I Casner et all 120031 ) on extended two-phase fluids and presented in Fig. 



[T]for wetting films and growing drops, the experimental conditions were to- 
tally different in the sense that instability was triggered by the total reflection 
of light within the deformation. Here, the beam incidence has been chosen 
to precisely avoid this mechanism, showing that an optically stretched drop 
can still become unstable above a radiation pressure threshold due to finite 
volume effects. We also investigated drop squeezing versus the amplitude of 
the radiation to Laplace pressure ratio, and for various drop radii. At large 
radiation pressure, the drop shape shifts from oblate to concave. As in the 
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stretching case, the squeezing significantly varies with the beam waist. With 
a further increase of the radiation pressure, the concavity reaches the solid 
boundary giving birth to a stable torus-like shape. Finally, by extending 
the developments on finite volume electrohydrodynamics to the optical do- 
main, our approach raised numerous questions on analogies and differences 
in drop deformations control by electromagnetic fields in general. Not all of 
these questions received definite answers. The optical jetting, for instance, 
could open new horizons in microdroplet dispensing. Consequently, even if 
the formation of stable cones and torus from radiation pressure effects de- 
serve experimental investigation, the present numerical study illustrates the 
opportunities offered by laser waves to actively manipulate droplets at the 
micrometer scale and reach dynamic and stable intriguing nonlinear drop 
morphologies. 
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